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A spherical cosmological model with a local void on scales of ~ 200 Mpc and with an 
inhomogeneous Hubble constant was proposed in recent two papers. This model explains 
consistently the observed properties of the cosmic bulk flow, the accelerating behavior of type 
la supernovae and the CMB dipole anisotropy without invoking a cosmological constant. As 
we are in a position deviated from the center in the model, the anisotropy of the Hubble 
constant appears owing to the directional difference of the distance from the observer to the 
boundary of the void region. It is found that the anisotropy is maximally about 6 % of 
the constant in the region of 200 ~ 500 Mpc from us. This inhomogeneity and anisotropy 
of the Hubble constant are not so large as to be inconsistent with the present observation. 
The detection of this anisotropy in the future will be useful to clarify the implication of the 
inhomogeneity of the Hubble constant. 



§1. Introduction 

A spherical cosmological model with a local void on scales of ~ 200 Mpc was 
proposed in our recent papers, 111 " ^ (cited as Paper 1 and Paper 2, respectively), to 
explain consistently the observed properties of the cosmic bulk flow, the accelerating 
behavior of type la supernovae (SNIa) and the CMB dipole anisotropy, without in- 
voking a cosmological constant. In this model the Hubble constant is inhomogeneous. 
Here we will derive the anisotropy of the constant associated with the inhomogeneity 
and discuss its consistency with the present observed values of the Hubble constant. 

The deviation of local expansion rates Hi, from the global Hubble constant Hq 
due to small-scale perturbations and large-scale structures of the universe has so far 
been studied by many people from various viewpoints. Here, let us briefly review the 
long history of these studies. First, Turner et al.Q) derived the probability distribu- 
tion function of Hq with given Hi, using numerical simulations, and Suto et al. ^ and 
Nakamura and SutoB* obtained it analytically in a simple model of inhomogeneities 
with a local low-density region. It was found as a result that the probability of a 
large deviation such as {Hi, — Hq)/Hi, > 0.4 is very small for inhomogeneities with 
the radii ~ 100/i~^ Mpc, where Hq = lOO/i km Mpc^^. In the linear analyses 
using power spectra in the CDM models, moreover, Shi and TurnerE* and Wang et 
al.@) showed that the cosmic and sample variances of {Hi, — Hq)/Hq are ~ 0.4% for 
the samples extending to ~ 100/i~^ Mpc. 

Second, the behavior of CMB dipole anisotropy was discussed to derive the 
constraints Jip on structures with inhomogeneous expansion rates bv^Nakao et a\..W , 
Tomita© -0)0), Shi et al.B, Shi and Turnerlzl', Humphreys et al.B), and Wang et 
al.i), where the anisotropy was assumed to be measured by an off-center observer 
in spherical models with an inner low-density region enclosed by an outer homoge- 
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neous region. When we consider only the redshift arising from the observer's peculiar 
motion relative to the CMB rest frame, the dipole component is found to be unac- 
ceptably large for the radius of the low-density region ~ 200/i"^ Mpc. The present 
author, however, showedllPUlP by solving exactly the null-geodesic equations from 
the off-center observer to the last-scattering surface, that in the spherical models, 
not only the redshift from the observer's relative motion, but also the gravitational 
redshift are effective, and their main terms are cancelled out, so that the total dipole 
component is small compared with the component only from the relative motion, 
and reduces to about 1/10 of the latter component, as long as he is near the center. 
Because of this special situation, spherical structures on the scale ~ 200/i~^ with 
Hq/Hi^ 0.8 are not constrained by the observed value of CMB dipole anisotropy. 

Third, the models with inhomogeneous expansion rates have been studied in 
connection with the spatial distribution and the [m. z] relation of type I supernoyae 
(SNIa) by Wu et al,a> Wu et al.,B Zehavi et alB and Tomita. i Zehavi et alM 
investigated the local inhomogeneity on the scale of ~ 60h~^ Mpc around the Local 
Group as a model of the Hubble bubble. However, it contradicts with the observational 
result of Giovanelli et al. which shows that the Hubble flow is uniform in the region 
within 100/i~^ Mpc. 

Independently of their model, the present author proposed another spherical 
model (the cosmological void model), which consists of the inner homogeneous region 
(on the scale of ~ 200/i^^ Mpc) and the outer homogeneous region with different 
Hubble constants (Paper 1 and Paper 2). It was introduced to explain the puzzling 
situation in the cosmic bulk flow on the scale of ~ 150/i~^ Mpc that was observed 
by Hudson et al.0 and Willick.© Their observational results are not comnatible 
with the other observations (Dale et al.,ll2P Giovanelli et al.,Ey) and Riess et al.Bil)) in 
homogeneous cosmological models, as was discussed by many people in the Workshop 
of Cosmic Flows (S. Courteau et al. 1999), but they may be compatible with them in 
our inhomogeneous model (Tomita©). This is because (a) in the inner homogeneous 
region there is no peculiar motion between comoving observers and clusters, (b) the 
off-center observers find a systematic motion of comoving clusters in the inner region 
relative to the global expansion, and (c) the total CMB dipole anisotropy measured 
by the off-center observers is small, and so the off-center observers in the inner region 
feel as if they were in the CMB rest frame, in spite of their motions relative to the 
global expansion. Thus our inhomogeneous model may explain the puzzling situation 
in the bulk flow, though the probability of its realization is very small. 

It was also shown that this model can explain the accelerating behavior of high- 
redshift supernovae, without a cosmological constant. It does not contradict with 
Giovanelli et al.'s result, because the diameter of the inner region is much larger 
than the size of their observed region WOh~^ Mpc. 

The local values of the Hubble constant have recently been measured by the HST 
Kev Project (Sakai et all^*), the High-z Supernova Search Team,© and Tutui et 
al.cB using the common calibrations due to Cepheid stars. The measurement in the 
first group using the multi-wavelength Tully-Fisher relation is limited to the region 
within lOO/i^^ Mpc and the obtained median value is 71 km Mpc~^. The unit 
km s~^ Mpc^^ is omitted in the following for simplicity. The measurements in the 
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second group were done in the region reaching 400/i^ Mpc, and the median value 
is 64. Thos e ir i the third group were performed using the CO Tully-Fisher relation 
(Sofue et al.E^EIP) in the region between 100 and 400/i~^ Mpc and the median value 
is 61. All of these values have errors of the order of ±10, but it seems that these 
locally measured values of the Hubble constant have the tendency to decrease with 
the increase of the distance. 

In view of the above studies we derive in this paper the effective Hubble constant 
(corresponding to the observed constant) in the above inhomogeneous models with a 
local void and study the anisotropy which will be observed by an off-center observer. 
The detection of this anisotropy will be useful to consider the implication of observed 
inhomogeneity of the Hubble constant. 



Various spherical inhomogeneous models (i.e. single-shell models, multi-shell 
models and models with a self-similar region) were considered in Paper 1. Here 
we treat only the simplest model consisting of the inner low-density homogeneous 
region (V^) and the outer high-density homogeneous region (V^^) connected with a 
single shell. The density parameters and Hubble constants in and V^^ are denoted 
as (I7g,i/g) and {Qq^.HJ^)^ respectively, where Hq = lOO/i' (/ = I or II). Here we 
assume Qq < and Hq > H^. The typical Hubble constant is Hq = 71 and 
i^ii = 57 (= 0.8 X Hi). The radius of the shell and the dist ance from the center C 
to our observer are assumed to be 200(/i^)~^ Mpc and 40(/i^)~^ Mpc, as in previous 



In this paper we consider the local behavior of light rays (reaching the observer) 
in the near region within 500 Mpc around O, and so the distance is expressed using 
the lowest-order terms of the expansion with respect to z{< 0.15). In constrast with 
the description in previous papers, the paths are expressed using polar coordinates 
{r,(j),6) with the origin O, as in Fig.l. 

The center C has the coordinates r = rc and (f) = tt. In these coordinates the 
line-element in and V^^ are expressed as 



where i = I or II and dQ"^ = dcjP' +sin^ cf) dO"^. The shell is given by = {(/)), which 
depends on time t', but its time dependence is neglected in the calculation of light 
paths, because {a} dr^^^ / dt^) / c (~ 10~^) is very small and the time difference for rays 
with (/) = and vr is also very small. 

In all rays reaching O are radial and straight, because of homogeneity in V\ 
and we have the relation 
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Fig. 1. A model with a single shell. 



between the origin and the boundary, where Hq = (dVa^)o- In the above equations 
we neglected the terms of 0{t^ — tj,)^- Then we have at the boundary 

c{tl-t\)=zlc/Hl (2-4) 

The direction of hght rays changes at the boundary in general, but the amplitude 
of their change is small and of the order of sinh r — r ?a ~ (-z^^)'^. which is neglected 
in the present approximation. Accordingly the light rays are regarded as straight, 
also in V^^, and so in V^^ we have the relation 



1 + z 



II 



1 + a"(t") 



or 



c(4^ - O = - ^P) c/Hl\ 



(2-5) 
(2-6) 



where we neglected similarly the terms of 0{t^^ — tY)"^- Because — Hq 
we can use H^^ in place of Hl^ in the present approximation, and so we obtain 

c{tY-t'') = {z''-zl')c/Hl\ (2-7) 

At the boundary the equality z\ = holds due to the junction condition that 
the lapses of times and t^^ are continuous (cf. Paper I). In the following we put 
z\{= zY), z^, z^^ as zi, z, z for simplicity. The dependence of the boundary is 
shown later. 

Accordingly the distance D between the observer and the source with redshift z 

is 

D {=c{tl-t^)) = z c/hI ioiz<zi{(f)) (2-8) 
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and 



D{= c{tl - t\) + c(tl^ - t")) = zi c/Hl 
+ {z - zi) c/H^ for z > zi ((/>). 



(2-9) 



At the boundary we have 



D^{ct>)=z^{4>)c/Hl. (2-10) 
Here we define the effective Hubble constant Hq^ , which corresponds to the 



observed Hubble constant, by 
Then Hq^ = Hq for z < zi, and 



D = z c/Hf. 



Hi 



3. 

Hi' 



11^ 



for z > zi. If we eliminate z using Eqs. (2-9) and ( |2-10D , we obtain 

Hf = Hl' + {Hl-Hl')D,{c^)/D 



(2-11) 
(2-12) 

(2-13) 



for z > zi. 

Next the functional form of Di (0) is given. In terms of the distance Db (between 
C and a point B on the boundary) and the distance Do (between C and O), we obtain 



^i(</') 



dI - dI sin^ . 



1/2 



Do cos ( 



(2-14) 



from the geometrical analysis in Z\COB in Fig.l. U D < D^ — Do ov D > Di, + Do, 
we have D < ot > Di{4>) for all 4>. If Db + Do > D > Di, - Do, we have D = Di(</>) 
for (j) = (pi specified by 

D2 _ jj2 



ixi = cos (pi = 



2DoD 



(2-15) 



where D^ = {D^ - Diy^. 

The expressions of the effective Hubble constant are given as follows (cf. Fig. 2). 

(1) For D<Db- Do, 

Hf = Hl (2-16) 

(2) For Db + Do> D> Db- Do, 

Hf = Hl fov(l)>(Pi (2-17) 

and otherwise Hq^ is given by Eq. ( 2-13| ). 

(3) For D>Db + Do, is given by Eq. (|]l|) similarly. 

In the case of (2) and (3), Hq^ has the minimum and maximum at = and vr, 
respectively. 

Now let us consider the angular average of the effective Hubble constant {{H(f^)) 
for (2) and (3). It is defined by 



(2-18) 
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Fig. 2. A surface with D = const, in the case of Dc + Do > D > Dc — Do. 

where 

/= / ' Mm[D,Di{(j))]dcos(j)/ f ' dcoscj) (2-19) 

J<t>a J <t>a 

for the average interval [(j)a^'Pb\- As examples we consider the whole-sky average 
{(f)a = 0, i^ft = tt), the northern-sky average = 0, = 7r/2), and the southern-sky 
average ((/>a = 7r/2, (/>ft = tt). Here the north is taken to be in the direction of C ^ 
O. The integrals (/) corresponding to their averages are denoted as and 
/s(i), respectively, in which i = 2 and 3 for the above cases (2) and (3), respectively. 
The expressions of these integrals are given in Appendix A. 

The maximum and minimum values of Hq^ (for = tt and = 0) and the 
values {Hq^)[ for whole-, northern- and southern-sky averages (i = A, N and S) were 
calculated in a typical example of = AO(h^)-^ Mpc, = 5Do, Hi = 71 and 

= 57, which were assumed in our previous papers. Their behaviors are shown 
in Fig. 3. 

As a result it is found that 

(1) The whole-sky average {Hq^)\ is constant in D < — Do, and decreases 
gradually from 71 (for D = - Do) to 62 (for D ~ 500{h})-^ Mpc), 

(2) The difference between {Hf)j^ and {Hf)s is about (2.5 ~ 1.3) for D = (160 ~ 
500)(/i^)-i Mpc, 

(3) The difference between the maximum and minimum values of Hq^ is about 
(4.0 ~ 2.5) for D = (160 ~ 5Q0){h})-^ Mpc. 



Anisotropy of the Hubble Constant in a Cosmological Model 



7 



75 



70 

u 

o 

X 

65 



60 

100 200 300 400 500 
D [(h')'' Mpc] 



Fig. 3. A diagram of the effective Hubble constant Hq^ (km Mpc~^) and the distance D (Mpc). 
The Hnes show the inaxinmm, the southern-sky average, the whole-sky average, the northern-sky 
average, and the minimum in the order from top to bottom. 



These behaviors of Hq^ seem to be consistent with the dispersive but decreasing 
tendency of the observed Hubble constant. 

§3. Concluding remarks 

If there is a spherical inhomogeneity in the Hubble constant, the anisotopy in it 
and the bulk flow relative to the global expansion necessarily occur, as long as our 
observer is not in the center. If this bulk flow is the one found by Hudson et al. and 
Willick, the corresponding anisotropy should be detected in the same direction, by 
the observations of nearby SNIa and galaxies (through CO Tully-Fisher method). 
Then the north is in the direction of the cosmic bulk flow, that is, Z = 260dr 15°, b = 
-1 ± 12° (Hudson et al. 1999), or / = 266°, b = 19° with la error (Willick 1999). 
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In this paper we assumed the simplest cosmological void models, in which the 
Hubble constant changes abruptly at the boundary. If we assume smoother models 
with a self-similar intermediate region, the change in the effective Hubble constant 
also may be somewhat slower. 

With the present assumption, considering only the lowest-order terms with re- 
spect to z, the dependence on the density parameter was neglected. The nonlinear 
full treatment with respect to z was given in Paper 2 and the [distance - z] relation 
was found to be consistent with SNIa data. 

Here a few comments are added as for our void. It refers to the low-density 
region (V^) with respect to the total matter density. The galactic number densities 
in the two regions and V^^ depend on the different complicated histories of galaxy 
formation, and their difference at the boundary may be rather small compared to 
the difference of total densities. 



Appendix A 

Integrals for the averaging 



The integrals (/) for the whole-sky average (^a = 0; ^6 = tt) for the cases (2) 
and (3) are 

/a(2) = ] [2L»(/ii + 1) + Do{Ji - J2)l , 



where 



1/2 



Dr 



Dr, 



-— sm — 1, 



Dr 



(A-1) 
(A-2) 



and 



J2 = fJi 



D. 



Do 
2 



+ (f I sin 



Dr 



1/2 



- (Ml)', 



/a(3) = -Do 



DcV 



1/2 



+ 



Dc 
Do 



sm 



-1 A, 

Dr- 



(A-3) 



(A-4) 



The integrals / for the northern-sky average (0a = 0, 05 = 7r/2) and the southern- 
sky average (0a = 7r/2,0ft = vr) are expressed as lN(i) and Isii), respectively, with 
i = 2 and 3. Here the north is taken to be in the direction of C — > O. 

For Dc> D> Db- Do, 



/n(2) = ^Do 



J1-J2 
Is{2) = D. 



+ Dni, 



For Db-Do> D> D^, 



/n(2) = -DoJi, 



(A-5) 
(A-6) 

(A-7) 
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/s(2) = + 1) - ^DoJ2. 



(A-8) 



For D> Db + D, 



/n(3) = 7a(3) - ^Do, 
/s(3) = Ia(3) + ^Do. 



(A-IO) 



(A-9) 
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